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Abstract 

There are two kinds of Tsallis-probability distributions: heavy tail 
ones and compact support distributions. We show here, by appeal 
to functional analysis’ tools, that for lower bound Hamiltonians, the 
second variation’s analysis of the entropic functional guarantees that 
the heavy tail q-distribution constitute a maximum of Tsallis’ entropy. 
On the other hand, in the compact support instance, a case by case 
analysis is necessary in order to tackle the issue. 
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1 Introduction 


During more than 25 years, an important topic in statistical mechanics theory 
revolved around the notion of generalized q-statistics, pioneered by Tsallis 
[1]. It has been amply demonstrated that, in many occasions, the celebrated 
Boltzmann-Gibbs logarithmic entropy does not yield a correct description 
of the system under scrutiny [2j. Other entropic forms, called q-entropies, 
produce a much better performance [2]. One may cite a large number of such 
instances. For example, non-ergodic systems exhibiting a complex dynamics 
0 . 

The non-extensive statistical mechanics of Tsallis’ has been employed to 
fruitfully discuss phenomena in variegated fields. One may mention, for 
instance, high-energy physics 0 - 01 . spin-glasses [5], cold atoms in optical 
lattices [6], trapped ions [7], anomalous diffusion i, i, dusty plasmas ra. 
low-dimensional dissipative and conservative maps in dynamical systems m, 
[T^ . [12], turbulent flows [12], Levy flights [12], the QCD-based Nambu, Jona, 
Lasinio model of a many-body field theory [TB], etc. Notions related to q- 
statistical mechanics have been found useful not only in physics but also in 
chemistry, biology, mathematics, economics, and informatics im, lEi, m. 

In this work we revisit the subject by appeal, in a classical MaxEnt phase- 
space framework, to the second variation of functionals. We find that such 
analysis guarantees a maximum of the Tsallis’ entropy only in the case of 
the heavy tail distributions. Our present treatment makes it advisable, on 
a more general MaxEnt framework, to always look at the second functional 
variation. We begin our discussion by remembering the concept of second 
variation. 


2 Second variation of a functional 

The essential concept that we need here is that of increment h of a functional. 
Note that the general theory of Variational Calculus has been developed in 
a Banach Space (BS) [20]. Particularly important BS instantiations are, of 
course, Hilbert’s space and classical phase-space. 

The MaxEnt approach in Banach space requires a first variation that should 
vanish and a second one that ascertains the nature of the pertinent extremum. 
This second variation is not usually encountered in MaxEnt practice, since 
one believes that the entropy possesses a global maximum. This second 
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functional variation is the protagonist of the present endeavor. The approach 
is described in detail, for instance in the canonical book by Shilov [20] (for 
local minima). It is simply explained. 

One needs to evaluate the increment h of a functional F at the point y of 
the Banach space one is dealing with. One has 

F(y + h) - F(y) = S'F{y, h) + ^S^F(y, h^) + e{h) (2.1) 


where 


lim 

h^O 


|g(^)ll 

\\h\\^ 


0 . 


( 2 . 2 ) 


By dehnition, 6^F{y,h) is the hrst variation of F (if it is linear in h). 
6‘^F{y, h) is F’s second variation, quadratic in h. If y is an extremum of 
F then 


S^F{y,h) = 0 (2.3) 

and it is a local minimum if 

S^F{y,h)>C\\h\\^ C>0, (2.4) 

or a local maximum if 

6^F{y,h) <C\\h\\^ C<0 (2.5) 

where C* is a positive constant and ||h|| stands for the norm of h. In phase- 
space, the object of our present concerns, 

||h||2 = J dfi (2.6) 

M 

where M is the region of phase-space one is interested in and /i the associated 
measure-volume for the concomitant space. We start our consideration with 
reference to the orthodox instance. 


3 Motivation 

The following considerations should motivate the reader to seriously consider 
the importance of elementary notions of functional analysis in q-statistics. 
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3.1 q-Exponentials as linear functionals or distribu¬ 
tions 

A generalized function (or distribution) is a continuous functional defined on a 
space of test-functions [21]. A typical such test space is the so-called JC space 
of Schwartz’, of inhnitely differentiable functions with compact support. 

One can prove [21] that x", dehned by 

xf = x", for x > 0, 

= 0, for X <0. (3.1) 

is a distribution possessing single poles at integers a = —k with residues (at 
the pole) 

= ( 3 - 2 ) 

with k = 1, 2,...., n,.... |21j . 

A function is a particular instance of a distribution, called regular distribu¬ 
tion. A singular distribution is that which cannot represented as a func¬ 
tion. For example, Dirac’s delta is such a singular distribution. Tsallis’ 
q-exponentials e^, dehned as 

eq{x) = [1 -f (g - l)x]+"", (3.3) 

becomes clearly a distribution dehned via 


eq{x) = [1 (g - l)a;])_ '' = [1 (g - l)a;] i-? if 1 -F (g - l)x > 0 

eq{x) = 0 otherwise. (3.4) 

Computations involving e„ should fruitfully appeal to distribution theory. 


3.2 An instructive example for second variations [[20 

Generally, the minimal condition 

5^F{y,h)>C\\h\\\ (3.5) 

with G > 0 cannot be naively replaced by the weaker restriction 

5‘^F{y,h)>.Q (3.6) 
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Consider, for instance, minimizing the functional 


1 


Hv) 

= J y‘^{x)[x — y{x)] dx. 

(3.7) 

It is easily seen that 

0 



y{x) = 0, 

(3.8) 

gives a functional extremum 

for F{y), that is, for y given by fl3.8P 

one has 


6^F{y,h) = 0. 

(3.9) 

The second variation 




5‘^F{y,h) = J xh‘^{x) dx, (3.10) 

0 

is > 0 for any function h{x) ^ 0. Thus, one may naively assume that fl3.8p 
yields a minimum for Fiy). 

To disprove such an assertion it is enough, given e > 0, to consider as y{x) 
any non-negative function that is positive at a: = 0, does not exceed e — x for 
X < e, and vanishes for x > e. For example, let y{x) = e — x ioi x < e and 
y{x) = 0 for x > 0. Then, 

e 

r A 

F{y) = {e — x)'^{2x — e) dx =— — \ (3-11) 

0 

For a y{x) = 0, F = 0, but the functional does not possess a minimum there. 
Similar considerations apply regarding local maxima apply if one considers 
the restriction 

6^Fiy,h)<C\\h\\^ (3.12) 

with C < 0, that cannot be replaced by the weaker condition 

[6^F{y,h)<0. (3.13) 

Choose, for instance, the functional 

1 

F{y) = j y^{x)[y{x) - x] dx, (3.14) 

0 

and repeat the above analysis. 
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4 Applying second variation 

4.1 Boltzmann-Gibbs’ Statistics 


In the general case the prior information consists of N mean values corre¬ 
sponding to the observables < Ri >, i = However, the points 

we are about to make here emerge already at the simplest level of just one 
observable, the Hamiltonian H (canonical ensemble). We limit ourselves to 
this instance in this work. Additionally, we assume that H is lower bounded. 
The MaxEnt variational functional becomes 



(4.1) 


with P the probability density MaxEnt is designed to encounter. H is the 
hamiltonian whose mean value is called < U >. Finally, a and 7 will repre¬ 
sent Lagrange multipliers. We consider now E^’s increment. 


/ {P+h) In(P-l-h) dfi+a 

f{P + h)H d/i- <U > 

M 

m J 


+ 


7 I J (P + h) dfi-lj + J Pln(P) dfi-a\J PH d^i- < P > | - 

\M / M \m 


We can also write 
Fs{P + h) - Fs{P) = 


7 1 y P d/i - 1 I . 

M 

(-1 -ln(P) + aP + 7 )h 


(4.2) 


W 


M 


dfi + 0{h?). 

(4.3) 

From (14.3p we hnd the hrst variation with its associated Euler-Lagrange 
equation plus the second variation as well. One has 


— 1 — ln(P) -I- aH - 1-7 = 0 


(4.4) 
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( 4 . 5 ) 


- jd^, < C\\hW‘ 

M 

From fl4.4p one gathers that 



M 

with Z the system’s partition function and [3 proportional to the inverse 
temperature. Eq. 04.51) yields 

- j ^ d^i = -Z J dfi<-Z j d/i = -Z\\h\\^ (4.7) 

MM M 

Notice that we can pass from the second to the third integral because 
is always greater than unity. This is a trivial point here, but not so when 
we consider Tsallis’ statistics below. Looking at 04.7p we see that one can 
choose C = —Z. Remark that these are classical considerations. Problems 
with 04. 6 p at T = 0 are thus not surprising, on account of Thermodynamics’ 
third law. As a bonus, we discover here that the bound-constant C is the 
partition function itself. 

4.2 Tsallis’ Statistics 

Here there is no single way of computing mean values for the theory [2]. 
Several options are available, that are today considered equivalent for all 
practical purposes [21 [22], because there is a “dictionary” that univocally 
relates two given probability densities Pi — P 2 , obtained using two different 
mean-values’ choices [22]. We consider in this work the three more important 
such choices and restrict ourselves to quadratic Hamiltonians. We insist in 
stating that the q-exponential is dehned as (x = {3H > 0) 

eq{—x) = [1 — (1 — if (1 — q)x < 1; 

= 0 otherwise (Tsallis' cutoff), (4.8) 

that tends to the ordinary exponential as g ^ 1. One speaks of long tailed 
distributions for [1 — (1 — q)x\ > 0 for all x > 0 and compact-support ones 
whenever the Tsallis cutoff becomes operative for some x—values. 
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4.3 Orthodox linear choice 

One evaluates mean values in the customary fashion, linear in P, i.e., 
< R >= Jj^RPdfi. The concomitant Tsallis el functional is 


Fs{P) = - I Pnng{P) dfx + a | / PH djj- < 17 > | +7 | j P dfx - 1 j , 

(4.9) 


M 


\M 


\AI 


For the increment we have 


Fs{P+h)-Fs{P) = - [(P+hyiugiP+h) dfx+a 

f{P + h)H d/i- <U > 

M 

La/ J 


7 


{P + h) dfx — 1 


La/ 


+ y \ng{P) dfi-alj PH dfi- <U> 

M \A/ 


7 I J P dfx — 1 

M 

Eq. fl4.10p can be recast as (see Appendix B) 


( 4 . 10 ) 


Fs(P + h) - Fs(P) ^ 


M 


q 


l-q 


P"-^ + aP + 7 


h dfx- 




qP^-'^'— d^i + 0{h^), 


( 4 , 11 ) 


M 


Eq. fl4.1ip leads to the following equations: 


l-q 


P«-^ + aP + 7 = 0, 


(4.12) 


qP^-^H d^i<C\\h\[ 


M 


(4.13) 














Eq. fl4.12p is the Euler-Lagrange one while fl4.13p gives bounds originating 
from the second variation. Thus, fl4.12p entails (using the procedure given in 

1231) 


a = /3qZ^ 


p = 


7 = T - 
g - 1 




Z= [I + - q)H]—^ d/i. 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


M 


For Eq. fl4.13l) we have, 


W = - I qP'^-^h^ diJ = - I qZ^-'^[l + /3(1 - q)H]'^^h^ dp. (4.18) 


M 


M 


In order to obtain a bound, we need to hnd a constant C (independent, in 
particular, of id). Thus, one needs to make sure that the bracket 


[l + (3{l-q)H]>0. 


(4.19) 


This entails 


0 < g < 1. (4.20) 

In such a case, the integral without the bracket is smaller or equal than the 
integral with the bracket and we have 

W <-qZ^-'^Wh^W =-C\\h‘^\\. (4.21) 

The present arguments, based on fl4.19p . guarantee an entropic maximum 
only for long-tail (or heavy-tail) Tsallis distributions. For compact support 
distributions, the present arguments are inconclusive. The maximum may 
or may not exist. One should further investigate things on a case-by-case 
fashion. For instance, if g > 1, the bracket [1 -|- /3(1 — g)ii] might remain 
positive for low enough jS. We will discuss this possibility in a different 
Section below. 
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4.4 Curado-Tsalllis mean values 

A second alternative way of obtaining mean values has been advanced in Ref. 
[23], where the authors dehne 


<u>= J pm dfi, 

M 

SO that the MaxEnt Lagrangian Fs becomes 


(4.22) 


Fs{P) = - j Pnng{P) d^i + a ( J pm d^i- < 17 > 1 +7 ( y Pd/i - 1 I , 

M \m / \m 

(4.23) 

and for the functional increment one writes (see Appendix B) 


Fs{P+h)-Fs{P) 


/ {P+hy lnq{P+h) dfi+a 

f{p + h)m dll- <u> 

\J 

M 

kJ 

A/ 


7 


{P + h) dfi — 1 


La/ 


+ y P" Ing(P) d/i - a j y pm dll- < P > j - 

M \m / 


7 I I P dfi — 1 \ (4.24) 

M 

Expansion in order of up to is now demanded for A) (P+h)'^, B) \nq{P+h), 
and the product A) B). Keeping only terms of order h and h?, one deduces 
from (I4.24P that the linear term in h yields 


q 


q- 1 

while de h^-term generates 


[l-a(g-l)P]P''-' + 7 = 0, 


- q[l- a{q - l)H]P^-^h^ dfi < C'||//|p 


M 


Using (I4.25P produces here (with the procedure of [23]): 

a = —{3 


(4.25) 


(4.26) 


(4.27) 


10 











(4.28) 


7 = 




P = 


q-1 


= e,i-(3H)/Z, 


Z= [l-(]{l-q)H]—. dfi 


M 


while fl4.26p leads to 


(4.29) 

(4.30) 


W = -fq[l-a{q- l)H]P^-^h^ dfi = 

M 

= -J + (g - l)(3H]^^+^h‘^ = 

M 

=-Z^-iq J[l +(3{q-l)H]^h^ dfi. (4.31) 

M 

Again, in order to obtain a constant bound we need the bracket in the last 
line above to be positive. This entails, again, heavy-tail distributions, here 
implying 


Thus, 


g > 1. 


(4.32) 


C = qZ^-'^. 


(4.33) 


The obvious demand g > 0 is satished given (I4.32p . As in the preceding 
subsection, the present arguments guarantee an entropic maximum only for 
long-tail Tsallis distributions. For compact support ones, these arguments 
are inconclusive. The maximum may or may not exist. One should further 
investigate things on a case-by-case fashion. For instance, if g < 1, the 
bracket [1 -|- /3(g — 1)H] might remain positive for low enough [d. We will 
discuss this possibility in a future Section below. 


4.5 Tsallis-Mendes-Plastino (TMP) mean values 

Following Ref. [23] we tackle the relationships 


<U>= 


f P^H dfi 

M _ 

/F9 d^l 

M 


(4.34) 
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Our Lagrangian reads now 


Fs{P) = - J Pnn,{P) dfi + a\J pm dfi- J Pm^iKU > \ + 

M \m M 


7 / F d/i - 1 


(4.35) 


Thus (see Appendix B), 


Fs{P + h)- Fs{P) = - j {P + hy ln,(F + h) dp+ 

M 


a 


/(F + hypdfi- hp + hydfi<u> 

+ 7 

I (F + h) dp — 1 

k) kJ 

M M 


La/ j 


+ 


J F" ln,(F) dfi-alj Pm d/^ - J P^fi < U > 

M \m M 


This simplihes to 


7 1 y F d/i - 1 I . 

M 

Fs{P + h) - Fs{P) = 


(4.36) 


q-1 


+ a{H- <U >) 


qpQ i ^ t 


M 


- [l-a{q- 1){H- < U >)]qP‘^-^h^ dp + 0{hy. 


(4.37) 


M 


Now we gather that 


q - 1 


+ a{H- <U>) 


qP<i-^ + 7 = 0 


(4.38) 
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(4.39) 


- i /[l - afe - 1){H- < U >)]qP«-^h^ < C\\h\\\ 

M 

From fl4.38p one finds, with the usnal procednre (see [23]): 


" l-/3(l-g) < P>’ 

(4.40) 

gZ^-''- [1 -/3(1 - g) < P >] 

(4.41) 

„ [1 -/3(i - 

z 

(4.42) 

Z = J[1 - (3{l - q)H]^ djj. 

(4.43) 


M 


A bit of algebra produces, from the above relations. 


W = -J[l- a{q - 1){H- < U >)]qPi-^h‘^ d^l = 

M 

= 11 dq, (4.44) 

M 

and setting , 


W 



[1 + (g- l)/3P]^/(''-i) 
l + /3(q-l) <U > 


h^djj,. 


(4.45) 


M 

One sees that, as in the two preceding instances, the bracket [1 + (g — l)/5hf] 
must be positive (long-tails!) so as to find a constant bound. This entails 


g > 1, 


and one has 


W < -qZ^-'^ j 


M 


that is 


C = q 


l + /3(g-l) < P> 
1 


(4.46) 


<C\\h\\^ (4.47) 


l + /3(g-l) < P> 




(4.48) 
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and we obtain the reqnired bonnd for the entropy to become maximal. More¬ 
over, from 04.471) ■ we get, once again, the here rednndant condition g > 0. As 
in the preceding two subsections, here our arguments guarantee an entropic 
maximum only for long-tail Tsallis distributions. For compact support dis¬ 
tributions, the present arguments remain inconclusive. The maximum may 
or may not exist. One should further investigate things on a case-by-case 
fashion. The comment made below Eq. 04.33P above is also pertinent here. 


5 Example: the Harmonic Oscillator (HO) 

Consider the simple Hamiltonian (in phase space) H = + Q‘^. We will 

reconhrm the second variation functional restrictions encountered above in 
the concomitant three q-statistics cases. 

5.1 Linear constraint < U > 

We start by remembering fl4.1ip . that required g < 1. One has 

- j dfx = - j qZ^-'^[l + /3(1 - g)(p2 + dfi < 

M M 

-J qZ^-i[l+(3{l-q){P^ + Q^)]ph‘^ dfi<-qZ^-'i J dfi = 

M M 

since 

[l + /3(l-g)(p2 + g2)]^ >1. 

We reobtain the restriction on heavy-tail Tsallis distributions 

5.2 Curado-Tsallis non linear constraints 

We recall (4.35) and (4.36). We had here the restriction g > 1 and deal now 
with 


(5.1) 

( 6 . 2 ) 


- j qll-a(q-l)(P^+Q^)]P'>-^h^ dfi = - J Z"-’g|l+/l(i,-l)(F"+Q^)]|-‘/i=o!^< 

M M 
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(5.3) 


- j qZ'‘-‘K‘ dfi = 

M 

1 

Since g > 1 we have [1 + (3{q — 1){P^ + > 1. Heavy tails once again! 

5.3 TMP constraints 

Here we must go back to (4.46) and (4.49). The operative restriction is g > 1. 
Accordingly, 


j[l- a{q - 1)(P2 + Q2_ < u >)]qP‘^-^h^ d/i = 

M 


qZ^-1 

1 + P{q-1) <U> 


J[l+P{q-l){P^ + Q^)]l-^h^df^< 


M 

qZ^-1 

1 + P{q-1) <U> 


\\h\\^ < C\\h\\\ 


We need again to appeal to heavy tail distributions. 


(5.4) 


6 The compact support instance 

We now consider the case of compact support probabilistic distributions in 
the formulation of Curado-Tsallis (similar arguments can be made for the 
other two possibilities). In such a case we need to satisfy, for a maximum, 
the relation 


-Z2-'?gy'[l + /3(g-l)P]i^h2 d/i < C'||h||2 (6.1) 

M 

A maximum is not guaranteed if 0 < g < 1. Consider now in more detail such 
a q-interval. H must be bounded by above in all phase space. By choosing 
fd sufficiently small, the bracket above does take a minimum positive value 
A and then we have 


<C||h||2 


( 6 . 2 ) 
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Thus, selecting 


- = C (6.3) 

we conclude that entropy does exhibit a maximum. Another way of viewing 
this argument is to consider that H have a maximum value R. Consider 
g < 1. More specihcally, g = 0.5. Then our critical bracket reads 


[1 - 0.5/3i?] > 0, 

entailing, for jS = l/kT, 


(6.4) 


T > 


R 


(6.5) 


i.e., a minimum temperature in order to guarantee the desired maximal con¬ 
dition that concerns us here. One might wish to speculate that, for lower 
temperatures, since no entropic maximum is possible, equilibrium might not 
be reached. 


Let us, for instance, H be given by 

H = {P^ + Q^)n{P^ + Ql-P^ + Q2) (6.6) 

where R is the Heaviside’s step function. We need [1 -|- /9(g — 1)H] to be 
positive. 

Selecting g = | and {3 = 1/kT small enough we have: 

0<P^ + Q^<P^ + Ql<^, (6.7) 

and thus 

T > {pS + Ql)/2k, (6.8) 

A = + (6.9) 


7 Conclusions 

Our second variation bounds are given, for the three Tsallis’ cases, by the 
C-bounds (I4.2ip . (I4.33p . and (I4.48p . respectively. Note that, since Z vanishes 
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at T = 0, we can not guarantee there a finite bound C, as required by the 
second variation protocol. 

The three bounds yield exactly the same conclusion; entropic maxima are 
guaranteed only for long tail Tsallis distributions. For compact support 
distributions, the present arguments are inconclusive. The maximum may 
or may not exist. One should further investigate things on a case-by-case 
fashion, as we have done, for particular instances, in Section 6. 

Summing up, the three Tsallis treatments, that were proved to yield identical 
predictions for mean values in [22], still give the same results concerning the 
requirements for entropic maxima. 

It is almost trivial to show that, for lower bounded Hamiltonians, a quantum 
n—levels treatment yields identical conclusions (see Appendix A). 

Our present treatment makes it advisable, on a more general MaxEnt stand¬ 
point, to always look at the second functional variation. 
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AppendixA 


We consider a quantum system with n discrete levels of positive energies e*, 
probabilities pi, and increments hi. The probability-vector P and increment- 
vector h belong to Consider, for instance, the orthodox linear choice for 
mean values, i.e., one evaluates mean values in the customary fashion, linear 
in the probabilities. The concomitant functional is 


. 2=1 


Fs{P) = > ) + 7 ( I , (A.l) 

i=l \ i=l 

For the increment we have 

n 

Fs{P+h)-Fs{P) = - '^{pi+hiY \ng{pi+hi)+a 


i=l 


y^jPi + hi)ei- <U > 


i=l 


+ 


7 


'^{Pi + hi)-l 


2=1 


2 = 1 


. 2=1 


X^Pnng(Pi) - a ( X^PiC- <u > 

n ^ 


7 


, 2=1 


Eq. (IA.2P can be recast as 


Fs{P + h)-Fs{P) = Y^ 

2=1 



pj ^ + aei + -f 


hi 


(A.2) 


<h‘j 




2=1 


Eq. flA.Sp leads to the following equations: 


1-q 


P‘1 ^ + ae, + 7 ^ 0, 


2=1 


(A.3) 


(A.4) 

(A.5) 
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Eq. flA.4p is the Euler-Lagrange one while flA.Sp gives bounds originating 
from the second variation. Thus, flA.dp entails 


a = /SgZ^ 

(A.6) 

7 = 

(A.7) 

g - 1 

_[1 + /?(! 

Vi 2- 

(A.8) 

n 

(A.9) 


2 = 1 


For the bound given by flA.SD we have 


^ = = - q)Q]^hl (A.IO) 

i=l i=l 

To obtain a constant bound, independent of the e*, we must demand [1 + 
/3(1 — q)€i] > 0 (long tail distribution!). Accordingly, 


fE = - qZ^-'^ll + /3(1 - g)6j < -qZ^-<^ ^ = -qZ 


2-q 


2=1 


2 = 1 


(A.ll) 

From (lA.lip we see that C = —qZ^ ^ and g > 0. As g < 1 we have finally 
for g the bounds 0 < g < 1. 
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AppendixB 

Consider a functional of P called Fs{P), given by 

Fs(P) = - f Pnn,(P) d^ = ^+ [-P 

J q-1 J 1- c 

M M 


dfi. (B.l) 


Thus, we can write 


Fs{P + h) = - I {P + hyingiP + h) dfi = j^ + 


M 


M 


(P + h)^ 

1 - q 


d/i. (B.2) 


Accordingly, 


Fs{P + h) = -- + 

q-1 


pg + gp(q-i)/i+£(£^pg-2/^2 , ,, , , 

- - -^-d/i + Oh^ B.3 

1-q 


M 


20 









References 

[1] C. Tsallis, J. of Stat. Phys., 52 (1988) 479. 

[2] C. Tsallis, Introduction to Nonextensive Statistical Mechanics Ap¬ 
proaching a Complex World (Springer, NY, 2009). 

[3] A. Adare et al., Phys. Rev. D 83 (2011) 052004. 

[4] G. Wilk, Z. Wlodarczyk, Physica A 305 (2002) 227. 

[5] R. M. Pickup, R. Cywinski, C. Pappas, B. Farago, and P. Fouquet, Phys. 
Rev. Lett. 102 (2009) 097202. 

[6] E. Lutz and F. Renzoni, Nature Physics 9 (2013) 615. 

[7] R. G. DeVoe, Phys. Rev. Lett. 102 (2009) 063001. 

[8] Z. Huang, G. Su, A. El Kaabouchi, Q. A. Wang, and J. Ghen, J. Stat. 
Mech. L05001 (2010). 

[9] J. Prehl, G. Essex, and K. H. Hoffman, Entropy 14 (2012) 701. 

[10] B. Liu and J. Goree, Phys. Rev. Lett. 100 (2018) 055003. 

[11] O. Afsar and U. Tirnakli, EPL 101 (2013) 20003. 

[12] U. Tirnakli, G. Tsallis, and G. Beck, Phys. Rev. E 79 (2009) 056209. 

[13] G. Ruiz, T. Bountis, and G. Tsallis, Int. J. Bifurcation Ghaos 22 (2012) 
1250208. 

[14] G. Beck and S. Miah, Phys. Rev. E 87 (2013) 031002. 

011109. 

[15] G. Wilk, Z. Wlodarczyk, Phys. Rev. Lett. 84 (2000) 2770. 

[16] J. Rozynek, G. Wilk, J. of Physics G 36 (2009) 125108. 

[17] G. M. Gell-Mann and G. Tsallis, Nonextensive Entropyinterdisciplinary 
Applications (Oxford University Press, New York, 2004). 

[18] S. Abe, Astrophys. Space Sci. 305 (2006) 241. 


21 


[19] S. Picoli, R. S. Mendes, L. C. Malacarne, and R. P. B. Santos, Braz. J. 
Phys. 39 (2009) 468. 

[20] G. Y. Shilov: Mathematical Analysis (Perganion Press, NY, 1965). 

[21] I. M. Gel’fand and G. E. Shilov, Generalized Functions, Vols. 1-2 (Aca¬ 
demic Press, NY, 1964). 

[22] G. L. Ferri, S. Martinez, A Plastino, J. of Stat. Mech. P04009 (2005). 

[23] C. Tsallis, R. S. Mendes, A.R. Plastino, Physica A 261 (1998) 534. 

[24] E.M.F. Curado, C. Tsallis, J. Phys. A 24 (1991) L69; Corrigenda: 24 
(1991) 3187 and 25 (1992) 1019. 


22 


